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Exact integrability and algebraic Bethe ansatz of the small-polaron model with the open 
boundary condition are discussed in the framework of the quantum inverse scattering 
method (QISM). We employ a new approach where the fermionic /^-operator which consists 
of fermion operators is a key object. It satisfies the Yang-Baxter equation and the reflection 
equation with its corresponding JT-operator. Two kinds of 'super-transposition' for the 
fermion operators are defined and the dual reflection equation is obtained. These equations 
prove the integrability and the Bethe ansatz equation which agrees with the one obtained 
from the graded Yang-Baxter equation and the graded reflection equations. 

1 Introduction 

There have been reported various one-dimensional (ID) integrable models. In ID systems, 
some physical properties are beyond the perturbative interpretation [IJ and therefore in- 
tegrable models have attracted much interest. Among the integrable fermion models, the 
simplest but non-trivial one is the small-polaron model, a spinless fermion model with a 
hopping term and an interaction term. It describes the motion of an additional electron in 
a polar crystal || and exhibits some properties that differ from those of the Fermi liquid. 

The integrability and other properties for the small-polaron model with open boundary 
condition have been studied in many ways. By use of the coordinate Bethe ansatz and the 
string hypothesis the critical value of attractive boundary potential for a boundary bound 
state was shown. [[| In the quantum inverse scattering method (QISM), the integrability 
and the Bethe ansatz equation were derived from the graded Yang-Baxter equation and 
the graded reflection equations. 0, [| || 

In the QISM, the graded Yang-Baxter equation plays an important role for the fermion 
models. The corresponding spin models through the Jordan Wigner transformation are 
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often useful. Using the correspondence, we obtain the graded Yang-Baxter equation from 
the Yang-Baxter equation of the spin model. [|7|, || There, the graded Grassmann product 
and the super-trace are introduced. We know that the small-polaron model is related with 
the XX Z Heisenberg model. 

Recently we introduced a fermionic i?-operator. It clarifies the integrable structure of 
the fermion models with the periodic boundary condition. The fermionic i?-operator satis- 
fies the Yang-Baxter equation. The 'super-trace' of the fermion operator is introduced and 
the commutativity of the transfer operator proves the integrability. |], [Hj Furthermore 
this fermionic i?-operator enables us to evaluate thermodynamic quantities for the fermion 
system through the quantum transfer matrix (QTM) method. In the QTM method, parti- 
tion function is estimated in two dimensional space (the quantum space and the auxiliary 
space). The transfer operator exhibits the integrable structure and the QTM leads to the 
partition function. The former operates on the quantum space, on the other hand, the 
latter operates on the auxiliary space. Because the fermionic i?-operator has the same 
footing on both spaces, the QTM method can be applied. |LT] 



In this paper we extend the fermionic i?-operator approach to the open boundary 
condition. In §2, we briefly review the periodic boundary case of the small-polaron model 
with this i?-operator approach, where 'super-trace' (Str) is introduced. The Yang Baxter 
equation leads to the integrability and the Bethe ansatz equation. We also investigate the 
properties of the fermionic /^-operator and define two kinds of 'super-transposition' (st 
and st). In §3, we prove the integrability with the open boundary condition. We obtain 
the reflection equation and the dual reflection equation and prove the commutativity of 
the transfer operator. In §4, we evaluate the eigenvalue of the transfer operator and get 
the corresponding Bethe ansatz equation. This result agrees with the one obtained by the 
graded Yang Baxter approach. The last section is devoted to the concluding remarks. 

2 Small-polaron model with periodic boundary con- 
dition 

We summarize the fermionic i?-operator approach to the small-polaron model with the 
periodic boundary condition. The Hamiltonian of this system is 

N 

H = / (1) 

H Jj+i = -t&j+i ] + c i c i+i} + Vn j n j+1 , (2) 

where t represents the overlapping integral and V the electron-phonon coupling. The 
fermionic i?-operator, 

R a b(u) = a(u){-n a n b + n a n b }] + b(u){n a n b + n a n b } + c(u){c\c b + c\c a }, (3) 

. . sm(u + 2rj) sin it . . . A , 

<u)= \ U , b(u) = ——, c(u) = l, 4 
sin 2r] sin 2r] 
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consists of fermion operators where n a = c\c a and n a = 1 — n a . This i?-operator satisfies 
the Yang Baxter equation, 

R 12 (u - v)R 13 (u)R 23 (v) = R 23 {v)R l3 {u)R l2 {u - v). (5) 

As is the routine for the quantum inverse scattering method, from the equation (|5|) we 
have the global Yang Baxter equation, 

Rab{u - v)T a (u)T b (v) = T b (v)T a (u)R ab (u - v), (6) 

where the monodromy operator is expressed by 

T a (u) = R aN . . . R al (7) 
= A(u)n a + B(u)c a + C(u)c{ + D(u)n a . (8) 

Here the suffix 1, 2 ... TV mean the quantum space and the suffix a, b the auxiliary space. 
This equation (||) leads to the commutativity of the transfer operators t{u), 

[t(u),t(v)] = 0, (9) 
r(u) = Sir a T a (u)=A(u)-D(u). (10) 

Here 'super-trace' for the fermion operators (Str) is defined as 

Str a X a = a (0|X a |0) a - a (l|X a |l) a , (11) 

where |0) a and |l) a are the fermion Fock states satisfying |l) a = <4|0) a and |0) a = c a |l) a . 
This means that the system has a sufficient number of conserved operators and therefore 
is integrable. The logarithmic derivative of the transfer operator at u = gives the 
Hamiltonian defined by ([!]) and (|2]) with 

f = -^V, V = ^t. (12) 
sin 2rj sin 2rj 

The energy spectrum is obtained from the eigenvalue of the transfer operator. The eigen- 
state \4>) and the eigenvalue A(w) of the transfer operator are given by 

M 

r(u)|0=A(«)|0, |0 = I]S(«i)|O>, (13) 

3=1 

where |0) is a pseudovacuum state. There is a constraint on Uj(j = 1 . . . M) (the Bethe 
ansatz equation) so that (|T3"D should be satisfied, 

/ x n N M , x 

b i U i) J .J~ a \ U 3 ~ U i) 

3¥=l,3=l 
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In terms of the solution {uj} of (|T5D, the energy spectrum is given by 

j M 

E=^-\og A(u) | w=0 = N cot 2 V - T -■ \ v (16) 

au z — ^ sin Uj sm[Uj + lr\) 

The fermionic /^-operator has important properties which are useful further in the 
analysis. We define two kinds of 'super-transposition' for the fermion operators (st and st) 
so that (X^ a ) sta = X a for any fermion operator X, 

X a = An a + Bc a + Cc{ + Dn a , (17) 
= An a + Bc\ - Cc a + Dn a} (18) 
X?» = An a - Bel + Cc a + Dn a . (19) 

The unitarity and the crossing unitarity conditions are satisfied as follows, 
Unitarity 



sin 2 u 



Rn(u)R 12 (-u) = p(u), p{u) = 1 - - ; 2~z ) (20) 

sin li] 

Crossing unitarity 

Ru(u)^R 12 (-u - 4 V f> = p(u), p(u) = 1 - sm2(M + 2V \ (21) 

sin 2r] 

We also have a property, 

Rn(u)^ = R 12 (uf 2 . (22) 

In the next section, we shall discuss the fermionic i?-operator approach for the open 
boundary condition. 



3 Exact integrability with open boundary condition 

To treat the open boundary condition, we consider the Yang-Baxter equation and the 
reflection equations. For the fermionic -R-operator given by @, we look for the fermionic 
.fT-operator which satisfies the reflection equation, 

R 12 {u - v)K 1 {u)R 12 {u + v)K 2 {v) = K 2 {v)R 12 {u + v)K 1 {u)R 12 {u - v). (23) 

We find that the fermionic i^-operator is given by 

, . sin(w — t) sin(u + t) 

K a (u) = y ——^n a + \ J n a . 24 

sin t sin t 
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We have a more general form of the fermionic K-operator, but we here restrict the K- 
operator to be even operator. Then the global reflection equation, 

R a b{u - v)7 a (u)R ab (u + v)7 b (v) = 7 b (v)R ab (u + v)7 a (u)R ab (u - v), (25) 

is satisfied where the double monodromy operator is 

7 a (u) = T a (u)K a (u)T-\-u). (26) 

We also show the dual reflection equation, 

R l2 (-u-)K?(u) 8tl R 12 (-u + - 4r])K+(vf* = K+ (vf* R 12 (-u + - A V )K+(u) s ^R 12 (-u-), 

(27) 

where R a b(u) = R a b(u) staStb and u± = u±v. This dual reflection equation ( P7| ) will be used 
for a proof of the commutativity of the transfer operator. The dual .ff -operator is found 
to be 

K+(u) = sin(w + 2r] - t + )n a + sin(w + 2r] + t + )n a . (28) 

Next we prove the commutativity of the transfer operators with different spectral pa- 
rameters. We define the transfer operator as 

t(u) = Str a {K+(u)7 a (u)}. (29) 

We prepare some properties of the 'super-trace' and the 'super-transposition' for the 
fermion operators. When Y is an even fermion operator and X operates trivially except 
on a space, we have 

{X a Y ab } st « =Y^Xf\ (30) 

When X and Y are even fermion operators, we have 

Str ab {X ab Y ab } = Str ab {X%;Y^}. (31) 

When X and Y are even fermion operators and X operates trivially except on a and b 
spaces, we have 

{A ab Y ab l - Y ab A ab , {6Z) 
StY ab {X ab Y ab } = Sti ab {Y ab X ab }. (33) 

The relation t(u)t(v) = t(v)t(u) is to be proved. We begin with 

t{u)t(v) (34) 

= StT a {K+(u)7 a (u)} Str b {K b + (v)7 b (v)} (35) 

= Str a {K+(u) st «7 a (u) st «} Str b {K+(v)7 b (v)} (36) 

= Str ab [K+(u)^7 a (u) s ^K b + (v)7 b (v)} (37) 
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Because X-operator is even, 7 a (u) sta commutes with K b {v). From the crossing unitarity 
(PI), the identity operator can be inserted as 

= pV 1 Stx ab [K+(u) Bt "K+(v)R ab {-u + - Ar,)^ R ab (u + ) s ^7 a (u)^7 b (v)} (38) 
= pV 1 Str ab i{K+(uf«R ab (-u + - 4 V )K b + (vfTn^)Rab(u + )7 b (v)}^] (39) 
= P+- 1 Str ab [{K+(u) st °R ab (-u + - 4rj)K+(v)^} st ^{7 a (u)R ab (u + )7 b (v)}) (40) 



The above equalities are due to the properties fl3T|) and ( |3T]) respectively, where u± 
u±v, p± = p(u±), p± = p{u±). Then from the unitarity (f2~0|) and property fl32|), we get 



P+-VI 1 Str af) [{K+(«) st ^ ab (-^ + - 4?7)K+( V )^} st ^ 

i2«6(-u-)i^(ti-){a' a («)i2 a6 (« + )T 6 (i;)}] (41) 
pVV: 1 Str afe [{ J R ab (- M _) st ^" t "K+( M ) st " J R ab (- M+ - 4p)if>) s " t *} st ^ 

(42) 



Here the reflection equation ( p3|) and the dual reflection equation (|27|) can be applied, 

= P+-VI 1 Str a6 [{if+(t;)*»i2 a6 (-tt + - A V )Kt(u)^R ab (-u^ st ^} stbAa 

{7 b (v)R ab (u + )7 a (u)R ab (u-)}] (43) 

= P+-V: 1 StT ab [R ab (-u^){K b + (vf»R ab (-u + - 4?7)K+( M )^} st ^ 

{T b («) R ab (u+)7 a (u) }R ab (u-)) (44) 



For this expression, we do the other way around. The property (|3!2j) is used in the above, 
and the property (|33|) and then the use of the unitarity (|2~0"D give 

= pV 1 StT ab [{K+(vf"R ab (-u + - 4 V )K+(u) st «} st ^{7 b (v)R ab (u + )7 a (u)}] (45) 
= t{v)t(u). (46) 

Thus, the commutativity of transfer operator is proved. 

We relate the transfer operator t(u) to the Hamiltonian of this system. The derivative 
of the transfer operator is considered. We have K a (0) — 1, 7 a (0) — 1 and 

£(0) = Str a K+(0) = -2cos2?7sint + . (47) 

Using these relations, we obtain 

N-l 

f (0) = Str Q K+'(0) + Str{K+(0)R' aN (0)P aN } + Str K+{0) £ ^ +li (0)P i+lj 

iV-1 

+ Str a K+(0)K[(0) + Str a Jf+(0) £ P i+ yi^ + y (0) + Str o {^+(0)P oiV < JV (0)} (48) 
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Here ' means the derivative with respect to u. The Hamiltonian density Hjj+x, (0) with 
(|12|), is expressed as Hjj + i = i^- +lj -(0) Pj+ij = Pj + ijR'j + ij(0). Thus the Hamiltonian is 
given by 

# = ^(0)-V(0)-Str a iC(0)) (49) 

N-l 

= H ii+i ~ ~ cott(2ni - 1) + - cott + (2n N - 1) (50) 
2 2 

3=1 

4 Algebraic Bethe ansatz 

In this section, the eigenvalue and the eigenstate of the transfer operator are discussed. 
The double monodromy operator (|26|) is expressed by 

7 a {u) = A(u)n a + T>(u)c a + Q{u)c\ + V(u)n a . (51) 

Because the double-mono dromy operator is an even operator, 25 (m) should have the prop- 
erty of the creation operator. We construct the eigenvector algebraically. We introduce 
the eigenvector of the transfer operator by 

M 

|$>=n^-)|0>. (52) 

3=1 



The double-monodromy operator ( pT|) is also expressed with the monodromy operators (|8|), 

7 a (u) = T a (u)K a (u)T-\-u), (53) 
T a (u) = A(u)n a + B(u)c a + C{u)c\ + D(u)n a , (54) 
T-\u) = A{u)n a + B(u)c* + C{u)c\ + D(u)n a . (55) 

The above relations give A{u) . . . D(u) in terms of A(u) . . . D(u) and A(—u) . . . D(—u). We 
have the relations between A(u) . . . D(u) and A(—u) . . . D(—u) obtained from the global 
Yang Baxter equation @ with v = —u 

T^\-u)R ab (2u)T a (u) = T a (u)R ah (2u)T b \-u). (56) 

It is convenient to use D(u), 

f>(u) = sin(2w + 2r])V(u) - sin 2rjA(u) (57) 

instead of D{u). The eigenvalue of T){u) acting on the pseudovacuum |0) takes a simpler 
form than that of D(u), 

A( u )\0) = S ^±^ a (u) m \0), (58) 
V(u)\0) = _ ^2,sin(, + 2,-t) ^ r|Q) ^ (5Q) 
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The transfer operator is rewritten as 

t(u) = StvK+(u)7 a (u) (60) 
= sin(w + 2t? + t + )A(u) - sin(w + 2t? - t + )V(u) (61) 
sin(2^ + 477) sin(n - t+) + sin^ + 2g + t+) - 



sin(2« + 2r]) sin(2it + 2r/) 

The global reflection equation gives the relations among A(u) . . .D(u). For instance, 
we have 

= (63) 
ai6 2 S(w)yi(t>) = 6 lC2 e(t;)D(M) + c x b % <Z(v)T>(u) + 6ia 2 yi(u)S(u), (64) 
cia2$(ti)D(v) + 6i6 2 ©(w)3(v) + Cl c 2 A(«)S(v) = aia a !B(v):D(ii) + aic 2 A(v)3(M), (65) 
aic 2 S(w)D(v) + aia 2 A(w)B(?;) = cic 2 S(v)2?(m) + b 1 b 2 'B(v)A(u) + cia 2 yi(t;)S(w), (66) 

where a% = a(u — v), a 2 = o(w + tj), . . . and so on. Using i)(it) instead of D(u), we obtain 

i>(«)2(t,) = ^ + ^ + + _ sin 2?? sin2(n + 2,) 

sin(M — u) sin(u + rj + 277) sin(it — v) sin 2(i> + 77) 

sin 2?? sin 2w sin2(-u + 2r?) _,. . „. . ,„_ N 

H -, — r — K — l Ar'B(u)A(v), 67) 

sm(u + v + 2r))sm2(v + r)) 

„. . sin(w — u — 2r?) sin(w + v) _ , , „, . sin27? . ~ . . 

sm('U — v) sm(u + v + It]) sm{u + v + 277) sin 2(tj + 77) 

sin 2r? sin 2v _ , . „ , . . . 

+ — 7 ^ - -^{u)A(v). (68) 

sin(w — v) sin 2(t> + 77) 

These relations lead to the eigenvalue of the transfer operator, 

sin(2u + 477) sin(u — t + ) sin(u + t) 2N t— r sin(?j — Uj — 277) sin(7x + Uj) 



cinl9^/ 4- 9n I Gin f J- X 



sin(2w + 277) sin t fj^ sin(u — w^) sin(w + ^ + 2?y) 

sin(« + 277 + t + ) sin 2m sin(u + 2rj — t) ^ . ^ 27V -r-r sin(u + + 477) sin(w — uj + 2rj) 
sin (2m + 277) sin/; -j-Jr sin(M — Uj) sin(w + itj + 277) 

(69) 

and the Bethe ansatz equations (1 = 1,... M), 

Ja(li/)1 2Ar Sin(77; + 277 + t + ) Sm(ui + 277 — t) -pr Sin(tti + Tij + 47/) sin(M; — TXj + 27/) 

\ fe(itz) J sin(tt; — £+) sin(w; + £) sin(«j — Uj — 2rj) sin(it; + Uj) 

(70) 

From the relation the energy spectrum is 



£?= -(JV- 1) cot 2t/- -cott + -cott + + — 2V ^m2-q 

2 2 sin 477 ^ sintijSin(Mj + 277) 
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5 Concluding Remarks 



In this paper, we have extended the fermionic i?-operator approach for the small polaron 
model with the open boundary condition. The Yang-Baxter equation and the reflection 
equations are satisfied among the fermionic i?-operators and the fermionic iT-operators 
which consist of fermion operators. The properties of the fermionic -R-operators and the 
'super-trace' and the 'super-transposition' for the fermion operators play an important 
role. This approach gives a simpler and clearer description than the usual graded Yang- 
Baxter approach. The integrability and the Bethe ansatz equation are obtained in a unified 
manner. 

We discuss future problems related to this work. First, we solve the Bethe ansatz equa- 
tions numerically to get the physical quantities. Second, we further extend the fermionic 
i?-operator approach to the boundary impurity problems. Last but not least, we develop 
the quantum transfer matrix method for open boundary condition. We note that this 
problem remains to be open. 
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APPENDIX A: Comparison with the spin system 

The small polaron model corresponds to the XX Z Heisenberg model through the Jor- 
dan Wigner transformation. For the periodic boundary condition (p.b.c), the fermion 
system and the spin system are not exactly the same. The energy spectrum and its Bethe 
ansatz equation for the fermion system are given by ( JO^ ) and ([16]) while those for the spin 
system are given by 

i M 

E = ± log A(u)\ u=0 = Ncot2 V -J2 -■ \ y (A.l) 

au t— ^ sin Uj sm{Uj + lr\) 

a(ui) \ N TT ( a{uj-Uj) \ 



We see that the right hand side of fllS)) and ( |A.2|) are different in signs. It reflects the 



difference between p.b.c. of fermion system and that of spin system caused by non-locality 
of the Jordan Wigner transformation. 

On the other hand, the small polaron model with open boundary condition is equivalent 
to the XXZ Heisenberg model with boundary magnetic field and their energy spectrums 
are exactly the same. 
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